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Cosmological issues of the gravitino production by the decay of a heavy scalar field X are exam- 
ined, assuming that the damped coherent oscillation of the scalar once dominates the energy of the 
universe. The coupling of the scalar field to a gravitino pair is estimated both in spontaneous and 
explicit supersymmetry breaking scenarios, with the result that it is proportional to the vacuum 
expectation value of the scalar field in general. Cosmological constraints depend on whether the 
gravitino is stable or not, and we study each case separately. For the unstable gravitino with M3/2 ~ 
lOOGeV-lOTeV, we obtain not only the upper bound, but also the lower bound on the reheating 
temperature after the X decay, in order to retain the success of the big-bang nucleosynthesis. It is 
also shown that it severely constrains the decay rate into the gravitino pair. For the stable grav- 
itino, similar but less stringent bounds are obtained to escape the overclosure by the gravitinos 
produced at the X decay. The requirement that the free-streaming effect of such gravitinos should 
not suppress the cosmic structures at small scales eliminates some regions in the parameter space, 
but still leaves a new window of the gravitino warm dark matter. Implications of these results to 
inflation models are discussed. In particular, it is shown that modular inflation will face serious 
cosmological difficulty when the gravitino is unstable, whereas it can escape the constraints for the 
stable gravitino. A similar argument offers a solution to the cosmological moduli problem, in which 
the moduli is relatively heavy while the gravitino is light. 
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INTRODUCTION 



It is quite plausible that the universe once experienced 
the epoch where its energy density is dominated by the 
coherent oscillation of a scalar field. A typical example 
is the inflaton oscillation after the exit of the de Sitter 
expansion period in the slow roll inflation [l|. Another 
example is dilaton and moduli fields in superstring theo- 
ries. 

The oscillating field eventually decays, followed by the 
reheating of the universe. Then, the radiation dominated 
era (the hot big-bang universe) commences. At the same 
time, however, the decay of the oscillating field may pro- 
duce unwanted particles.^ Recently, it has been recog- 
nized that the gravitino production at the decay of a 
modulus field can cause cosmological disasters p, !1| . It 
turns out that the decay width into the gravitino pair 
is non-suppressed and is comparable to that into other 
particles which couple to the moduli field with Planck 
suppressed interaction. Thus, the branching ratio of the 
moduli decay into the gravitino pair is sizable: in fact it 
is typically at a percentage level, or even higher. If the 
damped coherent oscillation of the modulus field once 
dominates the energy density of the universe, the over- 
production of the gravitinos at its decay would cause 
serious problems to cosmology. The gravitino decay, if 
it is unstable, would spoil the success of the big-bang 
nucleosynthesis (BBN). In addition, these unstable grav- 
itinos decay into the lighter superparticles, resulting in 
the overabundance of the lightest superparticles (LSPs). 



^ In general, the moduli fields lead to cosmological difliiculties 
known as the moduli problem Q. 



These arguments push up the gravitino mass to a region 
disfavored from the naturalness problem associated with 
the weak scale Q. 

One should be aware that the problem of the gravitino 
overproduction may apply to other heavy scalar fields. 
In particular, the case of inflaton is important. After 
the epoch of de Sitter expansion driven by the vacuum 
energy, the inflaton field falls down to its true minimum 
and starts the damped coherent oscillation around the 
true minimum. Eventually it decays to reheat the uni- 
verse. In many inflation models, the mass of the infla- 
ton, or more precisely the mass of the oscillating field 
after the slow-roll inflation, is much larger than the weak 
scale. Thus, in the low energy supersymmetry, the decay 
into a gravitino pair is likely to be kinematically allowed. 
If the branching ratio is not negligible, then the grav- 
itino production at the inflaton decay will cause serious 
cosmological problems. 

In this paper, we would like to investigate the decay 
of a heavy scalar field into a gravitino pair in a general 
ground and consider its cosmological implications. In 
Sec. II, we shall first discuss the partial decay rate of 
the scalar decay into the gravitino pair. The decay am- 
plitude is proportional to the vacuum expectation value 
(VEV) of the F auxiliary component of the scalar field. 
We will estimate the VEV of the F and thus the decay 
width into the gravitino pair in a very general setting. 
Then, in Sec. Ill, we will discuss cosmological problems 
caused by the gravitino problem. We will carry out the 
study in both unstable and stable gravitino cases, and 
identify the region of the parameter space which survives 
various cosmological constraints. A particular attention 
is paid to the case where the scalar field shares the prop- 
erties with the moduli fields. We will show that the case 
faces a serious cosmological difficulty when the gravitino 



is unstable, but can survive the constraints for a light and 
stable gravitino. We will also point out a new window of 
the gravitino warm dark matter with the mass range of 
10 MeV to 1 GeV. In Sec. IV, we will draw our conclu- 
sions and also discuss implications of our results to the 
inflation model building. In Appendix, we will explain 
the details on the estimation of the abundance for the 
long-lived superparticles. 

While preparing the paper, we received a preprint fa|, 
which dealt with a similar subject. Our results agree 
with Ref. 5], where overlap. 



II. HEAVY SCALAR DECAY INTO 
GRAVITINOS 

Let us begin by discussing the decay of a heavy scalar 
field X into a pair of gravitinos. To avoid unnecessary 
complication, we consider the case where the mass of X 
is much heavier than the gravitino mass, Mx ^ -^3/2- 
Also we assume that X is a singlet under the standard 
model gauge group. These assumptions are relevant for 
later use. 

The scalar field X may decay into a pair of gravitinos: 



X -^ ?A3/2 + "03/2 



(1) 



The decay is induced through the interaction in the su- 
pergravity. As recently calculated in Refs. H,|j{, its par- 
tial decay width is given by 



To proceed, we now present a general expression for 
the VEV of the X's auxiliary field. ^ In the absence of 
D-terms, the scalar potential of the A/" = 1 supergravity 
is given 



^SUGRA = e« [g^G^'^Gj - 3) , 



(4) 



where the subscript i (i) represents the derivative with 

respect to a scalar field 4>^ {(j) ) and G^^ is the inverse of 
the Kahler metric. The first derivative of the potential is 
obtained as 

= GiVsUGRA 

-fe^ (Gi-jG^'^Gi + GjG^'^iGj, + GjG^'^Gj,,) 

= GiVsUGRA 

Noticing that the F auxiliary field is expressed, in terms 
of G and its derivatives, as 



F' = -G^^e^/^Gj, 



(6) 



we find 



"FJ 



l^^suGRA = ~e'^/2gyF^"-e«/2(l + V)G,jF 



-G^fkF'F'' , 



(7) 



3/2 



"3/2 



Mi 



2887r Ml, 



(2) 



where Mp ~ 2.4 x 10^® GeV is the reduced Planck scale. 
The coupling constant d^/2 is defined by the relation 



J3/2' 



^"3/2 
M^ 



^((G,^)-e«/^G 



X 



(3) 



where G is the total Kahler potential, G = K + \n\W\'^, 
with K and W being the Kahler potential and the su- 
perpotential, respectively. The subscript X {X) denotes 
differentiation with respect to the X (X) field and (• • • ) 
stands for the VEV. Note that the right-handed side of 
the above equation is the (canonically normalized) F- 
auxiliary component of the X field. Therefore, ^3/2 pa- 
rameterizes the supersymmetry breaking felt by the X 
field. The parameter ^3/2 is expected to be order unity 
for a moduli field. If this is the case, the decay into the 
gravitino pair can be significant, in particular when the 
decays into other particles are also mediated by Planck 
suppressed interactions. Then, gravitinos produced at 
the scalar decay would lead to cosmological disasters. A 
more general case with less gravitino production will be 
discussed later. 



withV = GiG'^Gj- 3. 

Consider first the case where the supersymmetry is 
spontaneously broken within the framework of supergrav- 
ity. In this case, in addition to the X field, we need an- 
other field Z which dominantly breaks supersymmetry. 
At the vacuum, the VEV of the Kahler metric can be 
canonically normalized as (Gjj) = Jjj. With this basis, 
the stationary point condition reads 

/ 9 

— ( ^TT^SUGRA 

= -(e''^'Gxx){F'')-{e^/'Gxz){F^) 

e«/2(l + V)) (F^) - {Gxzz){F^){F^) 

-{Gxxz){P'')(P'^) - {Gxxx){F'')(f'^) 
-{G^^j^){F^)(f'^). (8) 



Notice that 



(e«/2^ = M3/2 , {F^) o^ M3/2. 



(9) 



A similar analysis can be found in Ref. 6|. 



When the mass of the X field is much heavier than the 
gravitino mass, one finds 



3<=/2G 



XX 



/eG/2 


f 

Kxx ~ 
V 


+ 


Wxx 

w 


/e«/2 


Kxx - 
V 


Kl + 


Wxx^ 
W J 


(e^/2 


Wxx) + 


0(M3/2) 


Mx + 


0(M3/2) 







(10) 



In the above we have used {Kx + Wx/W) <C 1. 

Using Eqs. ©, ® and jTOJ), the VEV of the F- 
auxihary field can be approximately expressed as 



(F-) 



-^({Gxz)hh/^{F' 



+ {Gxzz){F^){F' 
Here, evaluation of (Gxz) can be done as follows: 

, Wxz WxWz' 



(11) 



{Gxz) 



K 



xz 



[kxz 



W 

Wxz 
W 



W W 



Kx 



Wz\ 
W I 



fl2) 



It is natural that Kx and Kxz have VEVs comparable 
to {X)? Furthermore, (W^z/M^) = 0(\). Assuming that 
{Wxz/W) is negligibly small, which is the case when the 
X field is separated from the supersymmetry breaking 
sector in the superpotential, we find 



{Gxz)^0{{X)). 



(13) 



On the other hand, we also expect that (Gxzz) — 
(Kxzz) — C'((^))- To summarize the evaluations given 
above, we conclude'* 



UF-) 



^3/2 

Mx 



\{X)\, 



(14) 



as long as there is no cancellation between the terms in 
the right-handed side of Eq. 111|) . This in turn implies 



3^3/2 



\{X)\ 



(15) 



in the Planck unit.^ 

Next, we would like to discuss the case where the su- 
persymmetric anti-de Sitter vacuum is uplifted to the 



^ This is the case when the Kahler potential takes the form of 

K = f{Z,Z)XX + ■■■ with (Z) = 0(1), for example. 
** For a moduli field, (X) = 0(1). The relation should simply read 

^ For example, in the case where W = W{X) + A^(Z + /3) 7| and 
the minimal Kahler potential, we find that (^3/2 = v^|(^)| in 
the Planck unit. 



Minkowski one by explicit SUSY breaking terms. This is 
indeed the case for the KKLT set-up [Q], where the ex- 
plicit SUSY breaking is originated from an anti-D3 brane. 
The scalar potential of the effective theory is then written 
in the form 



y = VsuGRA + Vbroaking ■ 



(16) 



where Vsugra is the scalar potential of the supergravity 
given in Eq. Q. Vbreaking is the explicit SUSY breaking 
term which is of the form 



breaking — V f {X , X) , 



(17) 



with some real function f{X,X). 

The stationary point condition of the scalar potential 
(|16|l with respect to X now reads 



/dV\ ,y/df\ I d ^^ 

\dx/ ^ ^\ax/ + \ax^^"°^^ 



/oi\ 



= 0. 
It follows from the above that 

^ ' Mx \dX 



(18) 



(19) 



By imposing that the vacuum energy vanishes at the min- 
imum, we can determine the value of V^: V — iM'^,^/ {f). 
Then we obtain 

To proceed further, let us see a set-up of KKLT- type |3]- 
The function / is likely of the form |^, lOj 

ln/ = -nln(r + T-CC), (21) 

where T is an over-all modulus field and C represents 
a (matter) field, with n being some constant of or- 
der unity. With this form, we obtain {d\nf/dT) = 
0{{T)-^) = 0(1) and {din f/dC) ~ 0((C)), which im- 
plies (dlnf/dX) ~ 0{{X)). Thus, we reproduce the 
same relation (|14|) even when supersymmetry is explic- 
itly broken.^ 

It is amusing to point out that the relation (|14|l can be 
understood in softly broken global SUSY. For example, 
let us consider the following Lagrangian: 



C 



x^x + {e^MiX^x + h.. 

-[W{X)]p + \i.c., 



e^e^M^x'^x 



D 

(22) 



" Estimation for an overall moduli case was given in Ref. Illl . 



with Ml and M2 the soft SUSY breaking parameters of 
order M3/2- Here X denotes a chiral supermultiplet and 
X'^ its hermitian conjugate, unUke the rest of the paper. 
We can ehminate the auxiliary field F-^ by using the 
equations of motion 



F 



X 



{W^+MiX) 



The resulting scalar potential is written 



V = \MiX + Wx\ 



M^XX. 



(23) 



(24) 



When X has only Planck suppressed interaction, dtot 
becomes of order unity, as far as the Born approxima- 
tion is valid. For example, when X is the heavy mod- 
uli field and couples to gauge supermultiplets through 
Planck suppressed interaction in the gauge kinetic func- 
tion, the decay rate of X is coniputed to be in the form 
Eq. (EHJ with rftot = C(l) H 11 • With Eq. ^, the 
reheating temperature is estimated as 



Tr = 5.9x lO^GeVfit, 



Mx 



IQiOGeV 



(29) 



The stationary condition leads 



= (^MiX + Wx) (Wxx) + Ml ( MiX + Wx 
~-Mx(F^) + M|(X). 

Thus, we obtain 



(F-) 



Mi 



M 



X 



\{X)\ 



Mi{X) 
(25) 



(26) 



which is similar to Eq. (|14ll . This argument implies our 
result is insensitive whether the origin of SUSY breaking 
is spontaneous or explicit. 



III. COSMOLOGY OF GRAVITINOS 

A. Gravitino abundance 

Let us discuss the cosmological implications of graviti- 
nos which are produced by the heavy scalar field X . We 
will consider the case in which X dominates the energy 
of the universe when it decays. This situation can be 
achieved if X obeys the coherent oscillation with a large 
initial amplitude and also it is long-lived. In this case, 
the decay of X reheats the universe. We define here the 
reheating temperature Tr by 



where we have used g^iTfi) = 200. It should be noted 
that the reheating temperature shouldjje Tr > 7 MeV 
to be consistent with the BBN theory 12], which means 
that the mass of X should be 



Mx > 1.5 X lO^GeV 



1 



(30) 



Furthermore, the branching ratio of the decay channel 
X -^ i/'3/2 + V'3/2 is 



B 



3/2 



_ ^3/2 

Tx 



"3/2 

36 dL 



(31) 



It is seen that B^/2 can be much smaller than unity if 
c^3/2 ^ 1 (and/or if dtot ^ !)• Indeed, B3/2 should be 
suppressed enough to avoid cosmological difficulties, as 
we will see below. 

Now we are at the position to estimate the gravitino 
abundance. At the reheating epoch, gravitinos are pro- 
duced by the decay process X -^ ^"3/2 +"03/2 as discussed 
in the previous section. The yield of the gravitinos pro- 
duced by this process is estimated as 



rX 



T, 



v^ ^3 

^3/2 - 2^3/2m 



(32) 



which is defined by the ratio between the gravitino and 
entropy densities. Moreover, gravitinos are produced by 
the thermal scatterings at the reheating. We denote this 
contribution by YJ/^- The total yield is then given by 



Tn 



90 



t^'^9*{Tr) 



V^xMp , 



(27) 



where gt{Tfj) is the number of the effective degrees of 
freedom at T = Tr, and the total decay rate of X is 
denoted by Tx- Notice that Tx is determined by the 
interactions of X to other light particles in the minimal 
supersymmetric standard model (MSSM), and so it is 
highly model dependent. To parameterize the ignorance 
of the strength of interaction, we introduce dtot to express 
the total decay rate 



X 



€^Ml 
Stt Ml 



(28) 



^^3/2 - y. 



X 

3/2 



y. 



TH 

3/2 



(33) 



Notice that Y3/2 remains constant as the universe ex- 
pands, as long as there is no additional entropy produc- 
tion, and we assume it in the present analysis. 

We should mention here that there are potential 
sources of gravitinos in addition to the above mentioned 
ones. First, gravitinos would be produced much before 
the decay of X such that the abundance remains sizable 
even after the dilution by the reheating. Second, grav- 
itinos would be produced by the decay X — > X -|- "03/2 
where X is the fermionic partner of the scalar field X. 
In order to open this decay channel, the large mass hier- 
archy between X and X is required (see, e.g., Ref. [13|). 




10^ 10^ 10'' 10^ 10^10^°10^^10^^10^^ 



Mx [GeV] 



is of order unity, K?^^ 



VBBN 

-^3/2 



10 



-15_ 



'3/2 



10-16 for M. 



3/2 



1 TeV and 



10"" for M- 



3/2 



10 TeV. The constraint 



disappears only when the gravitino mass is very heavy: 
-^^3/2 ^ 100 TeV. We will not consider such a heavy grav- 
itino, as the region M3/2 ~ 10^ — 10"^ TeV is excluded by 
the overclosure of the LSPs |J] and the allowed region 
with M3/2 ^ 10'^ TeV is disfavored as the solution to the 
naturalness problem inherent to the weak scale. 

The yield of the gravitinos produced by the thermal 
scatterings is estimated as [l^ |17| 



K 



TH 
3/2 



1.1 X 10" 



Tn 



IQiOGeV 



(34) 



which increases linearly as Tr increases.^ On the other 
hand, the yield of the gravitinos by the X decay is given 
by Eq. ^. 

Consider the case where the X field decays through 
Planck suppressed interaction, in which the total decay 
rate Vx is given as Eq. |(2H|) with dtot = 0(1). Then 



FIG. 1: Contour plot of ^3/2 in the plane of Mx and (^3/2 
when Tx is given by Eq. 1281 with dtot = 1- The solid lines 



correspond to Y'3/2 = 10 ^'' , 
to right, respectively. 



10" 



10" 



and 10"" from left 



^3^?-"^i»-"''»(to^ 



Fg^a - 8-9 X 10^6 B3/24 



M 



X 



lOiOQeV 



(35) 



In the following we will not consider these possible con- 
tributions. Finally, the decay of X generally produces 
superparticles, i.e. i?-parity odd particles, followed by 
cascade decay into gravitinos. In particular, this contri- 
bution may be important when gravitino is the LSP. We 
will come back to this point later. 

From now on, we will derive the cosmological con- 
straints on the gravitinos produced at the reheating by 
the X decay. Since the constraints strongly depend on 
whether gravitino is unstable or stable, we will discuss 
each case separately. 



B. Unstable gravitino 

We first consider the unstable gravitino, especially 
when the gravitino mass lies in M3/2 ^ 100 GeV -10 
TeV as suggested by the gravity mediation models of 
SUSY breaking. The heavy gravitino with 1/3/2 ^ 1 
TeV decays soon after the BBN epoch. The decay rate 
of the gravitino into the MSSM particles is estimated 
as r ~ 193/(3847r)Af|y2Afp, corresponding to the life- 
time r ~ 2.4 X 10'^sec(lTeV/M3/2)^/^ In this case, the 
yield ¥3/2 is bounded from above, i.e. ^3/2 < ^^3^2^, 
in order to keep the success of the standard scenario of 
BBN 14]. Otherwise, the decay products of the grav- 
itino would change the abundances of primordial light 
elements too much and consequently conflict with the 
observational data. The recent analysis 15] shows that, 
when the hadronic branching ratio of the gravitino decay 



with the total abundance Y3/2 — ^3^2 +^3'/? • ^"^ ^^^^^ case, 
¥3/2 is determined by Mx and ^3/2, and Y,^,^ > ^3'/? '^'-"^ 
Mx < 1.4 X 10^°{B3/2/10-^^) GeV. In Fig.[T]we show 
the constant contours of 1^3/2 in the plane of Mx and 



23/2- 



Here we have set dt, 



1. It is seen that the mass 



Mx (and hence Tr) is bounded from above by the BBN 
constraint ¥3/2 — ^3'/^ < ^3^2^- Moreover, in order 
to avoid the overproduction of the gravitinos by the X 
decays, the coupling constant ^3/2 should satisfy 



d3/2 < 2.0 X lO-'dJ 



or equivalently 



B3/2 < 1.1 X 10" 



M. 



lOiOGeV 



Mx 



lOiOGeV 



't^BBN' 
^3/2 

10-16 



'VBBN^ 
^3/2 

10-16 



,(36) 



.(37) 



This gives a stringent constraint on the decay into the 
gravitino pair, or equivalently, the i^-component VEV 
of the X field. In fact, when the total decay rate is 
controlled by the Planck suppressed interaction (dtot — 
0(1)), this constraint excludes the case of ^3/2 = 0{1) 
as far as the gravitino mass is M3/2 ^ 100 GeV-10 TeV. 
Thus, the cosmological moduli problem cannot be solved 



^ We have neglected here the contribution from the helicity 1/2 
components of the gravitino. 
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FIG. 2: Contour plot of 13/2 in the plane of (1^/2 Mx and Tr. 
The solid lines correspond to Y3/2 = lO"^*^, lO"^'^, 10"" and 
10~^^, from inside to outside, respectively. The dotted line is 



the lower bound on Tr from -B3/2 < 1 for d; 



I3/2 



10" 



by simply raising the moduli mass |^ Q . At the same 
time, the above constraint poses a severe restriction to 
inflation models. In particular, the modular inflation sce- 
nario where one of the moduli fields is used as the inflaton 
would face a serious difficulty. 

The bound H3t)|l also implies that the scalar decay sur- 
vives the BBN constraint if d^/2 is small and/or dtot is 
large. The former can be achieved by some cancellation 
among the contributions in Eq. (|ll|l . or by the VEV of 
the X field {X) smaller than the Planck scale. The latter 
can be achieved if the X field has stronger interaction to 
the MSSM particles than the Planck suppressed one. 

To illustrate such a more general case, it may be con- 
venient to use Tji and take it as a free parameter. In this 
case, from Eqs. Q and (PTjl the branching ratio B3/2 is 
written as 



-B.s/2 — 



"3/2 



90 



^/^ 2887r \TT'^g^{TR)) T^Mp 



Mi 



(38) 



where -B3/2 < 1 as it should be. Then, ^3/2 is given by 
Eq. 134|l . while Y^,^ can be written by using Eq. (|38|l as 



Y, 



X 

3/2 



90 



192^ y^-^g^iTR)) TrMp 



(39) 



Quite interestingly, Y^,^ is inversely proportional to T^, 
whereas Yj^ is proportional to Tp. This means that 
the BBN observation puts not only the upper bound but 
also the lower bound on the reheating temperature. Fur- 
ther, it is seen that the total yield Y3/2 = Yr^,^ -\- Y^j^ 



can be determined from two parameters, i.e. d^^ji^x 
and Tp. In Fig. [21 we draw the contour plot of ¥3/2 
in the plane of d'^/2Mx and Tp. Here we have fixed 
9*{Tr) = 200, for simplicity.® It is clearly seen that 
Y^S^ puts the upper bound on d^^/^Mx as well as 
Tfj. This bound can be found as follows: In general, 



^3/2 - i^3/2 



'3/2MX 



equality holds when Y^^ = Y^,^, i.e., Tp ~ 1.2rf: 
from Eqs. (|Sl|l and l|H^ . It is obtained that 

Y;^^^ ~ 2.6 X 10^22 GeV"^ d3/2Mx , 
and hence Y^^^ < Y^^^ gives the upper bound 



'3/2 



'3/2 



Mx < 4 X lO^GeV ■ 



%/2 



'vBBN' 
10-16 



Notice that we have assumed Mx ^ M 



3/2 



(40) 



(41) 



1 TeV in 



this case. Therefore, when ^3/2 — 0{1), this gives a sever 
bound on the mass of the scalar field X. 

We should note that the parameter space in Fig. [21 
may contradict i?3/2 < 1- Indeed, the viable reheating 
temperature is 



Tn > lO^^GeV 



J_ / d3/2Mx Y 

A llO^GeVJ 
"3/2 



(42) 



As an example, we also show in Fig. [3 this lower bound 
on Tp for d; 



3/2 



10 



-10 



It is thus found that the bound 



is insignificant as long as ^3/2 is sufficiently large 
Finally, when ¥3/2 



B3/2 
find 



Y^^^ , the branching ratio is 
1.8 X lQ-*Mx/Mp and then from Eq. ^ we 



^3/2^3/2 < 2.9 X 10-" 



'vBBN' 
10-16 



(43) 



This again tells that, in order to avoid the overproduction 
of the gravitinos, we have to require (i) very small ^3/2 
and/or (ii) very small -B3/2. The former one may be 
realized by a VEV of X much smaller than Mp, while the 
latter one may demand interactions of X to the MSSM 
particles, whose strength is much larger than 1/AIp, to 
increase dtot- 

Next, we would like to discuss an additional constraint 
concerning the abundance of the LSPs. The LSP is stable 
when the i?-parity is conserved. To be consistent with the 
current observation of the cosmic microwave background 
radiation [13I, the present LSP abundance should be 



^Lsph^ <^dmh^ = 0.1051° ml ^ 



(44) 



* This choice is unsuitable for Tr < 1 TeV, however, our final 
results do not change much. 
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FIG. 3: Present abundance of the wino LSP S7^/i when 
Tx is given by Eq. 1281 with dtot = 1- The thick and thin 
sohd hnes are for M^ — 100 and 300 GeV, respectively. The 
horizontal dot-dashed line is the present dark matter density 
f^dm/l^ ^ 0.105. 



where h ~ 0.73 Is'l is the present Hubble constant in the 
unit of lOOkm/sec/Mpc. 

The LSP is produced by the decays of the gravitinos 
discussed above. We then obtain the upper bound on the 
yield of the gravitinos as 

^3/2 < 3.8 X 10-12 (^^) (JIL^) ' , (45) 
^'^ V 0.105/ VlOOGeVy ' ^ ' 



with Mlsp being the LSP mass. We can see that the 
bound is much weaker than the constraint from the BBN 
(^3^2^)' ^^'^ ^* gives no additional bound. ^ In addition, 
the LSP is produced directly by the X decay [T^ |23, 
w\ . Here we consider the neutral wino W as the LSP 
to maximize the annihilation cross section leading to the 
most conservative bound. The details of estimating the 
present abundance of W are found in Appendix „ 

In Fig. we show i^^h^ in terms of Mx by using 



"w 
We find that ^^h? becomes constant for 



Eq. ^. 

Mjf > 4 X 10^ GeV for M, 



w 



100 GeV, while n^h 



becomes larger as Mx decreases. (See the discussion in 
Appendix J. Therefore, we obtain the lower bound Mx to 
avoid the overclosure by the wino LSP, and Alx ^ 2 x 10^ 
GeV for Mt^ = 100 GeV. The bound becomes more strin- 



This is the case, provided that we consider the gravitino with 



M- 



3/2 



lOOGeV-lOTeV. If the gravitino mass is heavier, this 



gives a stringent constraint as discussed in Refs. |^ 



gent for larger My,. For example, Mx ^ 10^ GeV for 
Mvr, = 300 GeV. On the other hand, when we take Tr as 
a free parameter, ^^h^ depends on Mx as well as Tr. 
As estimated in Appendix„ the smaller values of Mx and 
Tr are excluded. This restricts the cosmologically viable 
parameters of the field X. Finally, we should notice that 
the bounds become severer when the LSP is composed of 
other neutralino components. 



C. Stable gravitino 

We now turn to the case of stable gravitino. This is the 
case when the gravitino is the LSP with exact i?-parity 
conservation. In this situation, the present abundance of 
the gravitinos is bounded from above in order to avoid the 
overclosure of the universe 22, .23j . For this reason, let us 
estimate the density parameter of the gravitinos, which 
is given by the present energy density of the gravitino di- 
vided by the critical density pcr- The present abundance 
of the gravitinos produced by the thermal scatterings is 
given by [iq 



"3/2" 



0.21 



Tr 



/ M: 



3/2 



10i"GeV/ VlOOGeV 



(46) 



Note that ^J/^ is dominated by the contribution from the 
helicity 1/2 components of gravitino and it is inversely 
proportional to M3/2. On the other hand, for the grav- 
itinos from the X decay, we find from Eq. 132() that 



"3/2 — 



Pci/sq 



0.027 / M3/2 \ _ 
h^ VlGeVy I 10 



yX 

^3/2 



(47) 



where sq is the present entropy density and Pcr/so 
3.6 X 10^^ h^ GeV. Thus, the total abundance is ^^3/2 



"3/2 



3/2 
must require 



rjJJJ. To avoid the overclosure by gravitinos, we 



f^3/2/i' < f^d 



(48) 



Especially, when ^3/2 — f^dm, the gravitinos constitutes 
the dark matter of the universe. As seen from Eq. H47|l . 
the bound on 13^2 in the present case is much weaker 
than that from the BBN for unstable gravitinos. 

Furthermore, we have to take into account an addi- 
tional constraint. Indeed, the gravitinos produced by the 
X decay face a constraint from the cosmic structure for- 
mation. Since the momentum of the gravitino at the pro- 
duction {p ~ Afx/2) can be much larger than its mass, 
its free-streaming at the epoch of the matter-radiation 
equality may erase the small scale structures which are 
observed today. This warm dark matter constraint leads 
to the upper bound on the present velocity dispersion of 
the gravitino vq |24| . The power spectrum inferred from 
the Ly-a forest data together with the cosmic microwave 
background radiation and galaxy clustering constraints 



puts severe limits |25|, |26|. From Ref. [2g we find ap- 
proximately Vq '^ 4 X 10~*. On the other hand, the 
present velocity of gravitinos produced by the X decay 
is estimated as 



^0=0 



1 f9*e{To)Y L 4M2/2 TaMx 



2 \9*s{Tr) 



Ml M3/2TR 



, (49) 



where Tq = 2.35 x 10"^^ GeV is the present photon tem- 
perature and g^tsiTo) = 43/11. Using Eqs. ^^ and (|T7|l . 
we find that 



qX ^ 3^3/2 M3/2TR 
3/2 - 2 MxiPcr/so) 

Therefore, vq can be written as 

wo^l.2xlO"^B3/2, 



(50) 



(51) 



where we assume that Vl^,^h? — Vtdrah^ — 0.105. Thus, 
the warm dark matter constraint on the gravitinos is 
translated to the upper bound on the branching ratio 



-S3/2 < 3 X 10" 



(52) 



It should be noted that this warm dark matter constraint 
can become weaken, as the portion of the gravitino dark 
matter in the whole dark matter density becomes smaller. 
Ref. [2^] gives the upper bound to ^^12 ^ 0.12 fidm to 
eliminate the warm dark matter constraint. Here we shall 
use this as a representative bound. 

Another potential constraint comes from the BBN, 
since additional energy from the gravitinos by the X de- 
cay may increase the Hubble expansion rate at T '^ 1 
MeV too much, which results in the overproduction of 
^He. The bound, however, is rather weak: ^3/2 ^ 0.35. 

Let us first consider the case where the X decays 
through the Planck suppressed interaction and the to- 
tal decay rate of X is given by Eq. H28|l . In this case, we 
find from Eqs. (gnj and |jl3) that 



njf^h^ ~ 3.9 X lO^^dtot 



M 



X 



10"GeVy vo-iGevy 



Mc 



^V 



Vtf/^h^ ~ 2.5 X lO^^t 



M^ 



B 



3/2 



10- 



lOiiQeV 



M. 



3/2 



O.lGeV 



(53) 



(54) 



^'^ A similar discussion can be applied for the gravitinos produced 
by the thermal scatterings. Although the constraint from the 
structure formation puts the lower bound on the gravitino mass, 
the interesting mass region which will be discussed here is far 
above this bound. 
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FIG. 4: Upper bounds on Mx in terms of M3/2, when Tx 
is given by Eq. 11281 with dtot = 1- The dotted lines corre- 
spond to the upper bounds from ^3/2/1^ < fidm^^ = 0.105 
for -B3/2 ~ 10~^, 10~* and 10~^ from left to right, respec- 
tively. The solid line is the bound from ^^3/2 < f^dm for 
B3/2 = 3 X 10~*, and the gravitino dark matter becomes vi- 
able above this line. The dot-dashed line is the upper bound 
on Mx from il,f/2 < 0.12Sldm avoiding the warm dark matter 



constraint in addition to ^3/2 < fid 



when B 



3/2 



10" 



Notice that Tr is determined from Mx as shown in 
Eq. H29() . In Fig. 0] we show the contour lines of 
^3/2^^ = 0.105 in the plane of M3/2 and Mx by vary- 
ing i?3/2. It should be noted that 11 JH ex M~}^ while 
ri?,2 0^ -^^3/2 • Therefore, for the heavier gravitino mass 
region, 12^,2 ^ ^J/2 ^^^ the gravitinos produced by the 
X decay contribute significantly the present energy of 
the universe. It is clearly seen that one can escape the 

overclosure of the if the mass Mx (and hence Tr) is suf- 

2 
ficiently small. The upper bound on Mx scales as M^,r^ 

for ^1^/2 — ^3/2 with smaller gravitino masses, whereas 
it scales as M~,'^^ for 1^3/2 — ^^3/2 with heavier gravitino 
masses. When B3/2 ^ 3 x 10~^, the upper bound on Mx 
becomes more stringent due to the warm dark matter 
constraint (see Fig.^J. 

We should stress here that the cosmological moduli 
problem can be solved in the small gravitino mass re- 
gion. We expect for moduli fields that S3/2 ~ 10"^ from 
^3/2 = C(l) and dtot = 0(1). Even in this case, there 
indeed exists the parameter space avoiding cosmologi- 
cal difRculties of gravitinos, where the gravitino mass is 
M3/2 < 0.1 GeV and the moduh (X) mass is Mx ~ 10^- 
10^ GeV, corresponding to 10^ GeV >Tr>7 MeV. Here 
the lower bound on Mx comes from Eq. (|30|) , whereas the 



upper bound is due to the warm dark matter constraint. 
The required large hierarchy between M3 /2 and Mx can 
be realized in a class of models of moduli stabilization. 
(See, for instance, Ref. z7].) 

Here we would like to point out that the heavy scalar 
X decay into a pair of gravitinos can oifer an interesting 
and alternative window of dark matter. When i?3 /2 ^ 3 x 
10"'', the free-streaming effect of the produced gravitino 
is small so that the gravitino can constitute the dark 
matter of the universe, i.e. ^^/2 = ^dm- ^"^ In Fig. ^ we 
present the contour line of ^^/2^^ — 0.105 with -63/2 — 
3 X 10"^ by the solid line. In the region above this line, 
the gravitino from the X decay can become the viable 
dark matter as long as ^^12 — ^dm- 

We should mention that the properties of the grav- 
itino dark matter can be different for different choices 
of parameters. Namely, the gravitino becomes the warm 
dark matter for -B3/2 ^ 3 x 10"**, and it gets cooler and 
eventually becomes the cold dark matter as the branch- 
ing ratio decreases. Furthermore, the mixed scenario of 
cold and warm dark matter is possible. In particular, 
the gravitinos produced by X can compose the warm 
component while those from the thermal scatterings can 
compose the cold one. More precise observations on the 
small scale structures of the universe enable us to test 
these hypotheses. 

The small branching ratio -B3/2 ^ 10"* indicates from 
Eq. (|31|) that d^/2 ^ O.ldtot- Though for the moduli fields 
a naive expectation will be ^3/2 ^ d^ot, the suppression 
of one order of magnitude may also be possible in some 
moduli stabilization mechanism. If it is the case, the 
decays of the moduli field with mass of C'(10^)-O(10^^) 
GeV will yield the gravitino warm dark matter. 

So far, we have assumed that the total decay rate of 
the X field is given by Eq. H28I) with dtot not very for 
from unity. Now we take T x or T^ as a free parameter. 
Notice that we find from Eq. (|38|l that 



M 



X 



2887r 
. ^3/2 , 



90 



R3 T^ M^ 

-°3/2 ^i?. ^'^P 



(55) 



which leads to 



^3/2^ 



4.1 X 10" 



'^3/2-53/2 
10-6 



r-SL 



V lO^GeV 



M. 



3/2 



IGcV 



(56) 



The total abundance Vl 



3/2 



^ "3/2 



riJ5 is then deter- 



mined by the two parameters, Tr and ^3/2^3/21 when 
the gravitino mass Af3/2 is given. 



The gravitinos produced by the thermal scatterings can be cold 
dark matter of the universe when Mo in < 1 MeV as long as 
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FIG. 5: Upper bounds on Tr in terms of M3/2. The dot- 
ted lines correspond to the upper bounds from Q,^i2h'^ < 



f^dm^ = 0.105 for d3/2B3/2 = W~ , W~ , 10""' and 10"* 
from left to right, respectively. The solid line is the bound 
from f23/2 < f^dm for ^3/2^3/2 = 3 x 10"*, and the grav- 
itino dark matter becomes viable above this line shown as 
the shaded region. The dot-dashed line is the upper bound 
on Tfl from 0,f/2 < 0.12f2dm avoiding the warm dark matter 



'3/2 
constraint in addition to ^23/2 < ^d 



when 0(3/2-63/2 = 10 



Similar to the case with unstable gravitino, we can 
obtain the upper bound on d^i2Mx to avoid the overclo- 

sure by gravitinos. Since fi3/2 > Vl^}^ = 2JVt^,^Vl^^, 



3/2 
^f/?'^^ < fidm/i^ ^ 0.105 results in the bound 



Mx < 5.7 X lO^GeV ■ 



-'3/2 



(57) 



Further, if the branching ratio becomes larger than 
Eq. (|5^ . this upper bound becomes tighter due to the 
warm dark matter constraint. We can see the bound is 
much weaker than Eq. (|41|l in the case of unstable grav- 
itino. 



In Fig. |31 we show the contour plot of ^1^/2 



Q 



dm 

in the plane of ^^3/2 and Tn by varying ^3/2-83/2. No- 
tice that the upper bound on -B3/2 H52|l can be ap- 
plied even in the general case, and hence we expect 
c?3/2^3/2 ^ 3 X 10"* since 0(3/2 < 0(1). In Fig. we also 
show n3/2h'^ ~ 0.105 with ^3/2-83/2 = 3 x 10"*. More- 
over, the upper bound on Tn from the warm dark matter 
constraint is also shown when ^3/2-83/2 = 10~^. We find 
that the allowed region is similar to that in Fig. 0] but 
the dark matter window for the gravitinos produced by 
the X decay is enlarged. 

Further, we show in Fig.|Slthe allowed parameter set of 
the gravitino dark matter for 0(3 /2 = 1 by taking the grav- 
itino mass as M3/2 = 1 GeV and 0.1 GeV. For smaller 
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FIG. 6: Allowed parameters for the gravitino dark matter 
with M3/2 = 1 GeV (the thick sohd hne) and A/3/2 = 0.1 
GeV (the thin solid line) when 0(3/2 = 1. Contour lines of the 
branching ratio are shown by the dotted lines, which corre- 
spond to B3/2 = 1, 10"^ 10"^°, lO"^'^ and 10"^° from below, 
respectively. The upper bound on -B3/2 from the cosmic struc- 
ture formation is shown by the dot-dashed line. 



^3/2, the corresponding branching ratio becomes smaller 
as shown in Eq. (|38|l . It can be seen that the gravitino 
dark matter is possible in a wide range of parameter 
space, and an important implication is that X should 
have the suppressed branching ratio of the decay into a 
pair of gravitinos as shown in Eq. I|52|) . 

Finally, we should argue that the NLSP (next-to- 
LSP) decay into gravitino restricts the allowed parameter 
space. Here let us consider, as an example, the case when 
a (right-handed) scalar tau r is the NLSP. The lifetime 
of r ^ V3/2 + T is 



Tr ~ 5.9 X 10 sec 



(My^y ( 



Mr 



V 



ViGevy viooGevy 



(58) 



and T decays during the BBN period or later for M3/2 ^ 
10 MeV when we take, e.g., Mr = 100 GeV. Therefore, 
when the gravitino mass is sufficiently small, the r NLSP 
becomes cosmologically harmless. On the other hand, 
when the gravitino mass becomes larger, the decay prod- 
ucts would dissociate or overproduce the light elements 
synthesized at the BBN epoch and would conflict with 
the observations. To avoid this, the number of t at the 
decay time should be small enough. In Appendix „ we 
estimate the abundance of r, Yf in terms of Mx and T^. 
The BBN constraint on Yr {MrYr more precisely) 
strongly depends on how r decays into hadrons \iB . llTl 
l2q |. According to the hadronic branching ratio Bh in 



Ref. [2g| and also to the BBN constraints for Bh = I and 
10"'^ in Ref. jT^I, we find that two regions, Tr ^ 10^ sec 
and 10^ sec ^ t^; > 10'^ sec, are cosmologically viable. 
Here we have used Yr ^ 10"^'^ by assuming that the 
reheating temperature Tr is higher than the freeze-out 
temperature of the stau annihilation, rp(see the discus- 
sion in Appendix) . Although we cannot find the relevant 
constraint for 10'^ sec ^ r? ^ 10^ sec (i.e. Bh = 10"'^- 
1) in the literature, we consider such a window is also 
allowed from the rough estimate of the BBN constraints 
interpolating in the region Bh = 10^"^-1. In this analysis, 
therefore, we take Tr ^ 10^ sec as a representative bound, 
which leads to M3/2 ^ 4GeV for Mr. We should stress 
here that our final conclusions leave intact even when the 
bound on Tr becomes severer. Since the lifetime of r de- 
pends on Af~?, the upper bound on M3/2 can be enlarged 
by larger Mr. Finally, this analysis is true for Tr > Tp. 
On the other hand, for Tr < Tp, the bound on ^^3/2 be- 
comes stronger or weaker in the low or high Mx region, 
respectively. The bottom line is that the BBN constraint 
on the r NLSP decay can be escaped for M3/2 < 0{1) 
GeV at least when the reheating temperature is higher 
than about Tp. The detail analysis in other parameter 
space will be done in future publication. 

The upper bound on the gravitino mass of Af3/2 ^ 
0(1) GeV implies that the maximal reheating tempera- 
ture is Tfj ~ 5 X 10^ GeV. Such a high reheating tem- 
perature is possible only when -83/2 ^ 10^^ such that 



a 



3/2 



nJS. For larger ^3/2 the upper bound on Tr is 
suppressed due to the gravitino production from the X 
decay. Note that the upper bound on Tr is directly trans- 
lated into the upper bound on Mx through Eq. H29|) . For 
instance, Tr < 5 x 10^ GeV gives Mx < 9 x 10" GeV 
for dtot = 1- (Note that Mx > 1.5 x 10^ GeV to have 
Tr>7 MeV.) 



IV. CONCLUSIONS AND DISCUSSION 

In this paper, we have considered the cosmological im- 
plications to the decay of the general heavy scalar field 
into the gravitino pair. Here we would like to summarize 
what we have obtained in this analysis. 

As was shown in the previous works [3, 3 , the decay 
amplitude to the gravitino pair is proportional to the 
VEV of the auxiliary component (F-^) of the heavy field 
X. We have thus presented the estimate for the VEV 
of F^ in a general setting: we have considered the gen- 
eral coupling between X and the fields responsible for the 
spontaneous supersymmetry breaking, and also we have 
considered the case of the explicit supersymmetry break- 
ing as well. In both cases, we have obtained the same 
and simple estimate for this value when only a single X 
field participates. The result shows that generally the 
VEV of F^ is proportional to the VEV of the X field, 
and thus the partial decay rate into the gravitino pair is 
suppressed when the X's VEV is smaller than the Planck 
scale. 
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We have then considered the various constraints on 
the gravitino production at the decay of the heavy scalar 
field X. The relevant constraints are different whether 
the gravitino is stable or not, and so we have discussed 
the two cases separately. In the unstable gravitino case, 
the constraints we have considered are 

1. The BBN constraint on the decay of the graviti- 
nos, producing hadronic as well as electromagnetic 
activities. 

2. The constraint on the abundance of the LSPs pro- 
duced at the gravitino decay. 

3. The constraint on the abundance of the LSPs pro- 
duced directly at the X decay. 

It is well-known that the BBN constraint puts the up- 
per bound on the reheat temperature Tr in order to sup- 
press the gravitino yield produced by the thermal scat- 
terings [iJI- In the case at hand, the decay of the X 
field produces the gravitinos as well. Since its yield be- 
comes larger when we lower the reheat temperature, the 
BBN bound gives the lower bound on Tr. We have iden- 
tified the allowed range of Tr both in the case where 
the X field is moduli-like and in the more general case. 
We have confirmed that if the gravitino mass lies in 100 
GeV-10 TeV range, no allowed region exists when S" is a 
moduli-like field, namely the field whose the decay into 
gravitinos is not suppressed {d^/2 — 0{l)) and the decay 
into other particles is controlled by the Planck suppressed 
interaction (dtot = C'(l)). On the other hand, for a more 
general case, we have found that the viable region of the 
parameter space does exist, but is severely constrained, 
as was shown in Fig. |21 

We have seen that the second constraint is weaker than 
the first one, provided that the gravitino mass is in the 
range given above. On the contrary, the third constraint 
excludes the case of very low reheat temperature, as the 
annihilation processes of the neutralino LSPs are not very 
effective there and thus the LSP abundance exceeds the 
observational abundance of the dark matter in the uni- 
verse. 

For the stable gravitino case, we have discussed the 
following constraints: 

1. The constraint on the gravitino abundance to avoid 
the overclosure of the universe. 

2. The constraint from the warm dark matter, namely 
the free streaming of the gravitino produced by the 
heavy X decay is small enough for the gravitino to 
be a viable warm dark matter. 

3. The BBN constraint on the decay of the NLSP par- 
ticles into the gravitinos. 

The first constraint is similar to the first one in the un- 
stable gravitino case, but numerically in the stable grav- 
itino case it is less severe. This makes a wider region of 



the parameter space cosmologically viable. In particu- 
lar the constraint on the gravitino abundance allows the 
moduli-like scalar field when the gravitino mass is lighter 
than 1-100 MeV, depending on the mass of the X field. 

The second constraint given in the above list also gives 
a significant constraint. We have found that it puts the 
upper bound on the branching ratio as -B3/2 ^ 3 x 10^ 
when the gravitino warm dark matter is the dominant 
component of the dark matter. This constraint disap- 
pears if the gravitino dark matter contribution less than 
12% of the total dark matter density. 

The third constraint is also quite stringent, but rather 
involved. With the lack of the complete analysis all 
through the relevant parameter regions, we have argued 
to put the upper bound on the gravitino mass (or the 
lifetime of the NLSP) as roughly of order 1 GeV when 
the NLSP is the stau weighing 100 GeV. 

In the rest of the paper, we would like to discuss the 
implications of our results to the inflationary scenarios. 
Assuming that there is no entropy production at a later 
epoch, our consideration gives a stringent constraint on 
the reheating process right after the inflation. When the 
infiaton is one of the moduli fields ^3i 1^ uM i ^^'^ the 
oscillating field is in fact moduli-like, then our results 
severely constrain the allowed gravitino mass region. In 
xirticular for the unstable gravitino, the analysis in Ref. 
can apply, leaving only the very heavy gravitino, e.g. 
for the wino LSP M3/2 ?J lO'^ TeV, which is disfavored 
as the solution to the naturalness problem on the weak 
scale. On the other hand, the case of a light and sta- 
ble gravitino becomes cosmologically viable (see Fig. 5). 
However the low reheat temperature is required to sup- 
press the gravitino abundance produced through ther- 
mal scattering, which may be inconsistent with a class of 
modular inflation with inflaton mass around 10^*^ GeV. It 
is interesting to note that a new window of the gravitino 
warm dark matter opens up in which the gravitino with 
the mass around 100 MeV, produced by the heavy moduli 
decay, constitutes the warm dark matter. Furthermore 
the reheat temperature is of the order 10^ GeV, which is 
a natural range we expect with the modular inflaton mass 
given above. The price we have to pay to realize this fas- 
cinating case is a slight suppression of the parameter 6^3/2 
by one order of magnitude, which, we suspect, should be 
possible within the framework of modular inflation. 

As an important remark, we would like to emphasize 
that the above argument to the modular inflation can also 
apply to the cosmological moduli problem. In the unsta- 
ble gravitino case, the cosmological constraints are too 
strong to be escaped unless the gravitino mass is heav- 
ier than, say, 10"^ TeV as was discussed above. On the 
contrary, the cosmological moduli problem can be solved 
when the gravitino is light and stable, say M3/2 ^ 0.1 
GeV. Such light gravitinos can be realized in the models 
of the gauge-mediated supersymmetry breaking. Com- 
pared to the modular inflation where the X mass is rather 
high, it is anticipated that the moduli mass is not very 
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far from the electroweak scale. Thus somewhat different 
region in the parameter space may be favored. For in- 
stance, the X mass much below 10^° GeV is also fine in 
this case. Anyway, the solution requires the large hierar- 
chy between masses M^i2 and Mx , which can be realized 
in a class of models of the moduli stabilization [23| ■ 

Let us come back to the implications to the inflationary 
models. Our results given in Section 3 indicate that, to 
make the inflaton decay cosmologically viable in a wider 
range of the parameter space, a smaller branching ratio 
of the decay into the gravitino pair (a smaller ^3/2) and 
a larger total decay rate (a larger dtot) is favored. In 
a simple class of the chaotic inflation, the Lagrangian 
is invariant under a Z2 discrete symmetry, a reflection 
of the field variable X -^ —X. Thus at the minimum 
of the vacuum, Fx will vanish and thus the decay of the 
inflaton into the gravitino pair does not take place and so 
the model does not suffer from the gravitino production 
problem. There are other inflationary models in which 
the parameter d^/2 becomes much smaller than unity, by 
realizing the VEV of X in an intermediate scale lower 
than the Planck scale. Examples include a new inflation 
model and also a hybrid inflation. Whether these models 
are really viable or not require detailed case study. A 
consideration can be found in Ref. |^. 

Finally to make dtot large, one can construct an infla- 
tion model where the inflaton interacts to other parti- 
cles with renormalizable couplings. An example is given 
as follows. Suppose an inflaton X couples to a pair of 
right-handed neutrinos in the superpotential as Vt^ = 
{fN/2)XN''N'^, where Jn is a coupling constant and the 
Majorana mass of N'^ is given by M^ = /^{X). Now, we 
set {X) = lO^^ GeV. ^^ in this_case, we find B3/2 
from the decay rate in Eq, 
oi q, -, N" + N" is Fatc : 
becomes much larger in some parameter region 



10~^ 

H28|) . The inflaton decay rate 
: MfjMx/{32n)/{X}^, which 
For ex- 



9 



ample, B3/2 ~ lO-^^ for Mx = 10^" GeV and Mn = 10 
GeV by taking 1^3/2 = {X)/Mp. Due to this fact, the 
cosmological constraints become weaker. 
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FIG. 7: 



Contour plot of the present abundance of the wino 



The dotted lines are f2^/i — 



LSP n^h^ for A% 

10, 1, . . . , 10~^ from left to right, respectively. 

GeV, a-j^ft^ ==: 1.3x10"^ 



For Tr > 10 
j^,„ —^.^,^^^ . The solid line shows n^/i^ = 0.105, 
and the left-bottom region is excluded by the overclosure by 
W. The horizontal dot-dashed line represents the lower bound 
on Tr{> 7MeV). 



APPENDIX 

In this appendix, we briefly explain an estimation of 
the abundance of the lightest superparticle in the MSSM, 
which is denoted by A. In the text, we consider A as the 
neutral wino for unstable gravitinos, while A as the stau 
for stable gravitinos. 

The decay of X and also the thermal scatterings pro- 
duce A, and then the number density of A, ua, can be 
found by solving the following coupled equations (see, for 
example, Ref. [33J and the references therein); 



Ackno'wledginents 



px + ^Hpx = -Txpx , 


(AT) 


PR + -iHpR = +Txpx , 


(A.2) 


fiA + "iHuA = {(7v) [nl^ - 


-nl) + ^iP^Px,{k.i) 



We would like to thank T. Moroi and A. Yotsuyanagi 
for fruitful discussisons. The work was partially sup- 
ported by the grants-in-aid from the Ministry of Educa- 
tion, Science, Sports, and Culture of Japan, No. 16081202 
and No. 17340062. 



^^ Such a VEV of the X field may be reah zed in the supersymmetric 
model of the new and hybrid inflation l29l . 



where px and pr are the energy densities of the X field 
and radiations, respectively, and the dot denotes a time 
derivative. Mx and Fx are the mass and total decay rate 
of X . Ba denotes the effective number of A per a X de- 
cay, {av) are the annihilation cross section which is ther- 
mally averaged, and rieq is the equilibrium value of the 
number density oi A. In these equations, H is the Hubble 
expansion rate which is given by H"^ = {px+ Pr) / {^Mp). 
The cosmic temperature is found from p^. Here we have 
neglected the co-annihilation effects on n^i, and the con- 
tribution to PR from A since it can be negligible, and 



13 



10 1=11111111 iiiiiiii iiiiiiiii iiiiiiii mill iiiiiiiii iiiiiiiii iiiiiiii iiiiH 



_ 10 



10° k 



!T- - Hr."-1'1 



> 

a 10"^ - 



_ 10" 



DC 



10' 



10' 



II mill 



-12" 



10 



10^ 



-lO''' / / ;■ / ,' / 
i' i iiriiiJ iiiiiJ iiiiiiJ iiuiiiil iiMiiiil iiiliiJ I mill 



1 0"^ 1 0^ 1 0^ 1 0^ 1 0^ 1 0^1 0^ °1 0^ ^1 0^^ 0^ ^ 
Mx [GeV] 

FIG. 8: Contour plot of VV in the Mx-Tr plane for M? = 
100 GeV. The dotted lines correspond to Yf = lO'^-lO"^*^ 
from left to right. For Tr > 10 GeV, Y? ~ 1.2 x 10"". 
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also assumed the rapid thermalization of A just after the 
production [Hill. 

In solving these equations, we assume that the energy 
of the universe is dominated by X for H > Tx , and we 
take the initial condition such that the maximal temper- 
ature of the dilute plasma (for H > Tx) becomes suffi- 
ciently high to keep A in the thermal equilibrium. The 
decay rate Tx is parameterized by the reheating temper- 
ature Tji as 



X 



90 



Mp 



(A.4) 



The abundance of A is determined from two parame- 
ters Tfi and Ba/Mx, in addition to (av) and Ma- For 



simplicity, we shall take Ba = 0.5 in this analysis, and 
the results for other values of Ba can be found by the 
rescaling of Mx as long as Tr is considered as a free pa- 
rameter. Then, Tp and Mx control the abundance of 
A, i.e. the yield of A, Ya, after A decouples from the 
thermal bath. 

When Tr > Tp {Tp is the freeze-out temperature of 
A and Tp ~ Ma/20), Ya does not depend on Mx. This 
is because A is still in the equilibrium after the decay of 
X completes, and its abundance is evaluated from (av) 
and AIa as usual. On the other hand, when Tr ^ Tp, 
Ya strongly depends on Tp and Mx. When Mx is suf- 
ficiently small, A is produced so abundantly by the X 
decays due to the source term in Eq. (jA.3p that Ya is 
determined by its annihilation effect at i? ~ Tx- In this 
case, Ya is inversely proportional to Tp. On the other 
hand, when Mx is sufficiently large, the annihilation at 
H r^ Tx becomes insignificant and Ya is determined from 
the contributions from the production by thermal scat- 
terings in the dilute plasma and also from the X decay 
aX H ^ Tx- In this case, we find that Ya oc Tr. 

Now we present our numerical results. First, we con- 
sider the case when A is the neutral wino W which is the 
stable LSP. In this case, we use l21l 



{av) 



M^ 
w 



2-k{2 - xw)' 



(A.5) 



where gw is the weak gauge coupling and xy^ = 
M'^/M— . In Fig. [3 we show the contour plot of the 

present abundance i^y^h"^ = M^Y^h? / {pcr I sq) by tak- 
ing M^ = 100 GeV. For Tp > 10 GeV, n^h^ takes a 
value of 1.3 x lO'^. 

When A is the stau r which is the NLSP, we use 35! 



{av) 



inai 



Mi 



(A.6) 



In Fig. |S1 we show the contour plot of the yield Y:f by 
taking Mr = 100 GeV. For Tp > 10 GeV, F? takes a 
value of 1.9 x 10~ 
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